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ABSTRACT 



The extremely rapid burst of TeV photons from Mkn 421 (15 May 1996) can 
be reconciled with the standard properties of a relativistic 7-ray emitting jet 
(bulk Lorentz factor T ~ 10; size ~ 10 17 cm) if one assumes that the electrons 
are accelerated in conical shocks with both opening and viewing angles ~ 1/T. 
If the injection time and the cooling time are much less than the photon crossing 
time, an emission ring moves along the jet and leads to the appearance of a very 
rapid flare, in satisfactory agreement with the observations. 



Subject headings: BL Lacertae objects : individual (Mkn 421) — galaxies : jets 
- gamma rays : observations — radiation mechanisms : non thermal 
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1. INTRODUCTION 

Several blazars, all of them associated with radio loud AGN, show an intense and 



variable 7-ray emission ( |Hartman et al. 1992|) . 



Various authors ( [Maraschi. Ghisclliiii fc Cclotti 1992 ; Dermer, Schlickciser fc 



Mastichiadis 1992| ; [Sikora, Begelman fc Rees 1994| ; |Inoue fc Takahara 1996|) have shown 



that the observed luminosity - up to ~ 10 48 ergs s -1 or more - and variability - typical 
time scales ~ several hours - can be understood if the emission is due to inverse Compton 
scattering of relativistic electrons on softer photons inside a relativistic jet. These jets move 
with bulk Lorentz factor T ~ 5 — 10, at a distance from the nucleus of order 10 16 — 10 17 cm: 
they appear to be the inner part of the well known radio jets, with similar T but larger sizes 
(10 19 — 10 20 cm). Since the 7-ray emission occurs closer to the central objects, it provides 
direct information about the physical processes effective in the central region of AGN. 

In this connection, the discovery of large (factor of order 20) and very rapid (time scale 
around 20 minutes) variations of the TeV emission from Mkn 421 on 15 May 1996 flGaidos 



et al. 1996| ) represents an important clue to the physical conditions in the central region 
of AGN. Variability timescales of 20 minutes would entail a source size ~ 10 15 cm (with 
problems for photon-photon opacity) or T ~ 10 2 (with problems for the total energy and 
the low probability of observing an event with a visibility cone ~ 1/T). 

Our purpose in this Letter is to present a model for the occurrence of rapid variations 
which uses the standard scenario but does not require extreme values for the Lorentz factor 
T or the size. The variability is related to anisotropies in the comoving frame of the jet 
associated with the geometry of the shocks responsible for the particle acceleration. We 
shall show that strong spikes in the very high energy spectrum are a natural consequence of 
a nearly conical shock geometry with characteristic opening « 1/T. 
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2. RELEVANT TIME SCALES 

We shall assume that the jet has a bulk Lorentz factor r > 1, scale length z and 
radius r; the scale length is defined as the distance over which the bulk properties of the jet 
change substantially: let $ be one such property, then z = \ d ^ dz \ ■ The observer's line of 
sight makes an angle 6 ~ 1/T to the jet axis. 

The time scales relevant to the variability are the following ones, as measured by a distant 
observer: 



1. the geometrical time needed to the photons to cross the emission region either in 
length (longitudinal time At z « z/cT 2 ) or in width (radial time At r ~ r/cT); 

2. the radiative time At em over which the emitting electrons lose a large fraction of 
their energy; if 7 is the typical electron Lorentz factor in the comoving frame, the 
dot indicates the time derivative, and the prime refers to quantities measured in the 
comoving frame, then At em ~ (j/jy/T; 

3. the injection time Aij n over which the electrons are fed to the emission region; as will 
be shown in the following, the radiative time is typically shorter than the geometrical 
time, and the electrons must be accelerated in situ; there are arguments ( BlandforcT 



fc Levinson 1995|) which point to the electrons being not only accelerated, but also 
generated in situ. 



From the point of view of the observer, a variability episode appears convolved with Ai in , 
At em , At r , and is dominated by the longest of them. On the other hand, in a diverging jet 
the emission mechanisms become very rapidly ineffective over distances longer than z, so 
that At z is an upper bound to the observed variability At D b s 
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At obs « mm[At z , max(At in , At em , At r )] (1) 

There are various arguments which suggest that r rs z/r; if this is indeed the case, and 
At r « At z (as we shall assume in the following), At f, s is necessarily of the order of the 
longitudinal time At z . In this paper we show how a plausible geometric arrangement might 
release this constraint, otherwise unavoidable. 

The injection time scale Ati n is very likely the shortest of all, of the order of several 
minutes, depending on the mass one chooses for the central black hole. The remaining 
time scales At em and At z are related in an interesting way to the jet optical depth 
to photon-photon absorption. This is due to the near equality of the cross sections 
involved, and to the radiation density in the jet being the target for both the emission 
and the absorption processes. The relation can be expressed in closed form under certain 
assumptions: only the jet contributes to the radiation field, the radiation spectrum is flat 
in uF u , the scattering process is always in the Thomson regime. Then one has 

p _ La T At em _ 3?r T 3 31 e u x 

~ mcV At z ~ U 7 ' Tll ~ IQtiT^ ej [ ] 

The quantity t is the compactness of the jet radiation field in terms of the observed 
luminosity L, r 77 is the jet photon-photon optical depth, and e u and are the upper and 
lower energies of the photon spectrum, respectively, in units of mc 2 . By multiplying the 
second and third of Eqs. |2| we finally get 



At z 77 ~ in^ 7 r { } 

As long as the right hand side is <C 1, one can have at the same time r 77 < 1 and a high 
radiative efficiency At em At z . In certain BL Lacs, however, as for instance Mkn 421, 
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the Klein-Nishina regime is reached and e u /^jT « 1; furthermore, if one requires that 
the jet is not loaded with too many pairs ( phisellini fc Madau 1996| , |Sikora et al. 1997]) , 
t 77 has to be -C 1, and a certain degree of fine tuning becomes necessary. In general, 
however, At em < At z . The models which have been fitted to typical 7-ray blazars have by 
construction At z equal to the observed variability time scale, which is several hours; Salvati, 
Spada & Pacini (1996), for instance, in the case of 3C 279 find r ~ 10 and z ~ 10 17 cm. 



3. GEOMETRY AND RAPID VARIABILITY 

The acceleration of the emitting electrons must occur in situ, and strong shocks are 
the most likely location. The orientation of the shocks with respect to the jet axis is 
preferentially at angles ~ 1/T. This is seen in numerical simulations ( |Bowman 1994] , |G6mez 



|et al. 1997] ), and is due to arcsin 1/T being the causal angle traced by a perturbation which 
moves isotropically at speed c in the comoving frame. Also, when the flow is deflected 
by the external medium, angles around 1/T are preferred, since at those angles the ram 
pressure levels off at values close to the internal pressure. Finally, it is plausible that At in 
and At em are much smaller than At z . 

We can envisage a geometry of the emission region similar to Fig. [I]: the accelerating 
shocks are a series of opening and closing cones; a very thin layer of extra electrons (the 



perturbation) move down the jet with velocity /3c, (3 — y 1 — 1/T 2 ; they radiate immediately 
after being accelerated, so that the emission region is a ring at the intersection between the 
electron layer and the conical shocks; if the cone opening is precisely arcsin 1/T = arcos/3, 
and the line of sight coincides precisely with one generatrix of the cone, the photons emitted 
in the observer's direction pile up at the same observed time. 

Many of these restrictions can be released. One can assume that the cone opening 
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angle a and the viewing angle 9 are different from each other; then the pile up condition 
becomes 



(3C C n COSa fA\ 

9 = a ± arcos — — (4) 



cosa cos(6* — a) (3 

Real solutions exist only if a > arcos/3; when a > arcos/?, the pile up effect takes place 
along two lines of sight. Also, it is not essential that 9 coincides exactly with the pile up 
value. Fig. shows the result of viewing angles differing by various amounts from the 
optimal one. If x, y, and z are cartesian coordinates of the emission ring, with z along the 
jet axis, s = z tana is the ring radius, t is the time, and t b s the time at which a given 
wavefront reaches the observer, the flux F is given by 



F(t obs ) oc 



dX ^+(^-) 2 dt (5a) 



\ \dt obs J t \dt 



obs 



x = -r—p\tobs — t(l — Pcosd)], y = a/s 2 — x 2 , z = cf3t (5b) 

olllf/ 



The curves have been computed with 9 = e x arcos/3 (where e is a number < e < 1), 
a = arcos/?, T = 10, and At z ~ 9 hours. One should note that the vertical axis is logarithmic, 
and the height of the spikes is very often a factor of several; one should further note that 
blazars are observationally selected to have 9 < arcos/5, and the a posteriori probability 
of having 9 > e x arcos/5 is relatively large, P(9 > ex arcos/3 | 9 < arcos/?) = 1 — e 2 ; 
for instance, the sharpest peak shown in Fig. |2| should be observable in 19% of the cases. 
The peaks appear because the relevant geometry is not the one of the jet as a whole, 
but instead the one of the emission region, which -in the case of strong radiation losses- 
traces the surface of the shocks. The preferential orientation of this surface results in the 
radiation field being strongly anisotropic at any point of the emission region at the time of 
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the perturbation; the peak photons are distributed over an angle much narrower than 1/r, 
of the order of (1/r) x ^ At em / At z ; the relevant length is the radiative length, instead of 
the geometrical one, and one has At f> s « At em = (z/cT 2 )(At em /At z ). Both results can be 
described in terms of an "effective" Lorentz factor wTx J ' At z j At em , and explain why -at 
variance with Eq. 0- At b s <C At z « At r . 

Anisotropies in the laboratory frame on scales smaller than 1/r are equivalent to the 
loss of isotropy in the comoving frame, at variance with the standard scenario, and may 
have interesting consequences. While the peak photons, as we have found, are very closely 
aligned, the emitting electrons interact with the same peak photons at angles of the order 
of 1/r; this breaks the assumptions under which Eq. |3] was derived, and allows a situation 
where at the same time the radiative losses are very rapid and the pair production opacity 
is very low. Also, if the line of sight is not closely aligned with the shock surface, the 
observer will infer a photon density at the source lower than the one experienced by the 
electrons; this breaks the assumptions under which synchrotron self Compton models are 
usually computed, and allows larger ratios between the Compton and synchrotron regions 
of the spectrum. 

We stress finally that the sharp features carry only a small fraction of the total radiated 
energy; most of this energy resides in the long lived, low level tail, which by construction 
has a time scale of order At z . So the observed light curve will exhibit conspicuous sharp 
features only if the injection events are spaced at intervals > At z ; in the opposite case a 
given peak will be swamped by the superposition of the many tails of the preceding peaks. 

Because of the assumption of negligible radiative lifetime, the curves of Fig. |2] are the 
Green functions of the variability problem for the given set of parameters, and by convolving 
them with an arbitrary injection function one obtains the general solution. Fig. |3] is our 
best fit by this method to the TeV light curve of Mkn 421 on 15 May 1996 with the same 
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parameters of Fig. |2| and e = 0.9. It shows that it is not necessary to assume extreme jets 
to account for extreme variability; only, similar episodes should not be observed to repeat 
before times of order At z ~ several hours. We recall that Mkn 421 was observed to flare 
with a ~ 60 minute timescale also on 7 May 1996 flGaidos et al. 19"96p . With a different 
injection function we can easily fit this flare as well. Because of the longer duration of the 
episode, we find that the ratio of the fluxes after and before the flare should have been 
higher than on 15 May 1996, but there are no data to check this. 



4. CONCLUSIONS 

We have developed a model to explain extremely fast variations in the emission of 
relativistic jets. Under the following assumptions: 

1. the jet has typical parameters (T »s 10, z ~ 10 17 cm, At z ~ 9 hours); 

2. cold electrons are injected in the center and move outwards with bulk velocity c/3] 

3. electrons are accelerated in conical shocks with opening angle ~ 1/r; 

4. injection time and cooling time are <C At 2 ; 

5. the viewing angle is 9 ~ 1/r, 

it is possible to observe a very short bump <C At z , because the component of the peak 
photons' velocity along the jet is near to the electron velocity, and the photons pile up on 
the same wavefront. In particular, we have been able to fit with typical jet parameters the 
20 minute TeV flare of Mkn 421 on 15 May 1996. 

Deviations from the ideal geometry could be some curvature in the shocks, or a 
non-negligible lifetime of the electrons at low energies; these should result in a certain 
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amount of smoothing, which needs to be computed in detail. The extension of this model 
to the case of non-negligible particle lifetime is currently under investigation with the aim 
of comparing the expected behaviour of the emission at different frequencies. We may 
note straight away that in at least one case (PKS 2155-304) the burst has been seen to 
evolve with frequency-dependent time scales ([Urry et al. 1997| ), roughly compatible with 
synchrotron cooling in a uniform field. This is possible only if At r <C At z , r z/T (with 
far reaching implications on the jet dynamics), or in a non-standard geometry as the one 
proposed here. 

This work was partly supported by the Italian Space Agency through grants 
ASI-95-RS-120 and ARS-96-66. 
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Fig. 1. — The geometry assumed for the acceleration and emission region. See text for 
details. 

Fig. 2. — The light curves observed for a = arcos/5 and e = 6*/arcos/3 < 1. The values of e 
are indicated by the labels. 

Fig. 3. — A fit to the TeV flare of Markarian 421. The geometrical time scale of the 
underlying model is At z = 9 hours. 



